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Abstract. We present a complete calculation of the QCD and SUSY-QCD corrections to neutralino 
pair annihilation into bottom quark-antiquark pairs through exchange of a pseudoscalar Higgs 
boson, which is the dominant process in the cosmological A-funnel region of the mSUGRA model. 
We present numerical predictions for the annihilation cross section and discuss the influence of the 
correction terms on the cold dark matter relic density with respect to recent cosmological data. 

PACS. 12.38.Cy Summation of perturbation theory - 12.60. Jv Supersymmetric models - 95.30. Cq 
Elementary particle processes - 95.35.+d Dark matter (stellar, interstellar, galactic, and cosmo- 
logical) 
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1 Motivation 



Thanks to the recent WMAP mission and further 
cosmological observations, the matter and energy de- 
composition of our Universe is known today with un- 
precedented precision pQ. Direct evidence for the ex- 
istence of Cold Dark Matter (CDM) is accumulating 
[2J, and its relic density J?cdm can be constrained to 
rather narrow range [3] 



To evaluate the number density n of the relic par- 
ticle, one has to solve the Boltzmann equation 



0.094 < HcuMh 2 < 0.136 



(1) 



at 95% confidence level, h denotes the present Hubble 
expansion rate Hq in units of 100 km s -1 Mpc . 

Although the nature of dark matter remains still 
unknown, extensions of the Standard Model (SM) of 
particle physics provide interesting candidates for these 
Weakly Interacting Massive Particles (WIMPs). In Su- 
persymmetry (SUSY), a natural candidate is the Light- 
est Supersymmetric Particle (LSP), which is stable, if 
i?-parity is conserved. It is usually the lightest of the 
four neutralinos, denoted Xi or shortly \- 

The Minimal Supersymmetric SM (MSSM) is gov- 
erned by 124 a priori free soft SUSY-breaking param- 
eters, which are often restricted to five universal pa- 
rameters that are imposed at the unification scale and 
can be constrained using data from high-energy col- 
liders. As the neutralino relic density also depends on 
these parameters, its computation is another power- 
ful tool to put constraints on the parameter space and 
provide complementary information, in particular at 
high energies and masses that would otherwise not be 
accessible at colliders. 



= -3-ffn - (o-cfff) (n 2 



(2) 



with the Hubble rate H and the thermal equilibrium 
density n eq . The number density n is directly related 
to the relic density f2cDMh 2 = m x n/p c cx (crctfv)^ 1 , 
where m x is the LSP mass, p c = 3H 2 / (8itGn) is the 
critical density of our Universe, and Gn is the grav- 
itational constant [4]. The effective cross section a e s 
involves all (co) annihilation processes of the relic par- 
ticle x into SM particles, and (a e sv) signifies the ther- 
mal average of its non-relativistic expansion (u < c). 

Several public codes perform a calculation of the 
dark matter relic density within supersymmetric mod- 
els. The most developed and most popular ones are 
DarkSUSY [6] andmicrOMEGAs [7]. All relevant (^anni- 
hilation processes are implemented in these codes, but 
for most of them no (or at least not the complete) 
higher order corrections are included. However, due to 
the large magnitude of the strong coupling constant, 
QCD and SUSY-QCD corrections are bound to affect 
the annihilation cross section in a significant way. They 
may even be enhanced logarithmically by kinematics 
or in certain regions of the parameter space. 

It is the aim of this work to study these correc- 
tions for the neutralino-pair annihilation into a bot- 
tom quark-antiquark pair through exchange of a pseu- 
doscalar Higgs boson A° . This process dominates in 
the so-called A-funnel region of the mSUGRA parame- 
ter space at large tan /3, which is theoretically favoured 
by the unification of Yukawa couplings in Grand Unifi- 
cation Theories (GUTs) jSJ. Supposing a WIMP mass 



Cosmology 



Contributed Talk 






Fig. 1. The tree level diagram for the process XX ~~ * ~~ * 
bb, indicating the two possible fermionic orientations of the 
Majorana initial state. 



of 50 - 70 GeV, this process has been claimed to be 
compatible with the gamma-ray excess observed by the 
EGRET satellite in all sky directions [9]. However, the 
corresponding scenarios may lead to antiproton over- 
production, so that they would not be compatible with 
the observed antiproton flux [TP] . 



2 Analytical results 

As we focus here on the phenomenological aspects of 
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this work, the analytical calculation is only briefly sketcheci nn ^ 1 ^ at * on P rocess XX 
For a more detailed discussion we refer the reader to 
Ref. HI]. 

In the calculation of the leading order (LO) am- 
plitude, we have to make sure to anti-symmetrize the 
neutralino pair in the initial state properly by tak- 
ing into account both possible fermionic orientations 
shown on the right-hand side of Fig. [2] leading to a fac- 
tor \[2 in the amplitude with respect to the amplitude 
of one fixed orientation. Denoting by \fs the centre-of- 
momentum energy and (3 X = v/2 and (3b the neutralino 
and 6-quark velocities, the leading order cross section 
can then be written as 



Fig. 2. QCD (left), top-quark loop (right top), and SUSY- 
QCD (right centre and bottom) correction diagrams for the 

A -> bb. 



where Z\qcd, <4 to p and Z\susy relate to the QCD, 



o-lov = - 



1 (3 b 



Ncg 2 Tl n h\ bbi 



2 8ns Is - 



(3) 



It is proportional to the inverse of the flux factor sv, 
the integrated two-particle phase space s[3b/(8irs), the 
number of quark colours Nc = 3 and the squares of the 
weak coupling constant g, a neutralino mixing factor 



top-quark loop, and SUSY-QCD correction diagrams 
shown in Fig. [2] The same relation also holds for the 
Higgs boson width Ta, for which the same correction 
diagrams are involved. Using standard methods, we 
compute the QCD correction Z\qcd at 0(a s ), coming 
from the gluon exchange and real emission diagrams, 
which agrees with the known result for pseudoscalar 
Higgs boson decays in the on-shell scheme [12] . In the 
limit m 2 -C s, this correction develops a logarithmic 
mass singularity, which can be resummed to all orders 
using the renormalization group, i.e. replacing vrib with 
the running mass m&(s) in the Yukawa coupling YiAbb 
[13j . The remaining finite QCD corrections in the MS 
scheme are known up to G(p? s ), 



^QCD 



a s (s)C F 17 
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N 2j - tan 8 w N l3 ) (N M cos (3 - N 3i sin 0) 
3), (4) 



+ 
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(35.94- 1.36n/) 



(164.14- 25. 77n/ + 0.2597iJ), (7) 



the bottom quark mass through the Yukawa coupling 
h>Abb = —gmb taa/3/(2mw), and the Higgs boson prop- 
agator. The Nij in Eq. ^ are the entries of the matrix 
N that diagonalizes the neutralino mass matrix. The 
non-relativistic expansion of Eq. (j4|) is performed by 
expanding the squared centre-of-momentum energy s 
in powers of v 2 , 



(5) 



where n / denotes the number of flavours [14j . 

A separately gauge- independent correction at order 
0{a 2 s ) is induced by the top-quark loop diagram shown 
in Fig. [21 Its contribution [15] 
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(8) 



= 4m 2 1 + - + 0(v*) 



and is in agreement with the results given in Ref. [5]. 

The next-to-leading order (NLO) annihilation cross 
section can be written as 



can be important at small values of tan/3. However, it 
is largely suppressed in the Higgs funnel region con- 
sidered here, as shown in the centre panel of Fig. [3] 

In SUSY, the sbottom-gluino exchange diagrams 
shown in Fig. [2] give rise to additional corrections at 
0(a s ). The self-energy diagram leads to the mass renor- 
malization [8] 



Arrib 



a s (s)C F nig 



cnlo = clo (l + ^qcd + Aop + A 



(A b - n tan (3)I{m 2 ,m? ,m?) 
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(6) 



(9) 
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where /(mi, ma, m.3) = Co(0, 0, 0; mi, m%, ms) is the 
3-point function at zero external momentum. In this 
low-energy (LE) limit, and neglecting Ab with respect 
to the tan /3-enhanced /i, the vertex correction equals 
the mass renormalization [16] up to a factor 1 / tan 2 /3, 
so that the total SUSY correction becomes 



.(LE) _ a s (s)C F 1 + tan 2 /3 
SUSY " tt tan 2 /3 
x I(m~ bi ,m~ b2 ,m s ). 



m.g/ztan/3 



(10) 



As Ab may be of similar size as fj, tan /?, its contribution 
must also be resummed. This is effectively done by 
replacing [T7] 



lim Amb 



1 + lim u t 



fj.tan/3—fO 



Anib 



(11) 



Finally, our result agrees with the those in Refs. [TB] . 
and we implement the finite 0(m,b, s, 1/ tan 2 (3) remain- 
der as described in Ref. [T7] . 



3 Numerical discussion 

For our numerical study of the impact of the QCD, 
top-quark loop, and SUSY-QCD corrections discussed 
above, we place ourselves in a minimal supergravity 
(mSUGRA) scenario with five parameters mo, m 1 / 2 , 
Aq = 0, tan (3, and sgn(/i) at the grand unification 
(GUT) scale. The weak-scale MSSM parameters are 
then determined through renormalization group run- 
ning using SPheno j!9j . and the physical Higgs and 
SUSY masses with FeynHiggs [20 . For the SM inputs, 
i.e. masses and widths of the and electroweak gauge 
bosons and quarks, the coupling constants, the entries 
of the CKM-matrix, and the CP-violating phase, we 
refer the reader to Ref. [21] . 

We first study the effect of the corrections on the 
annihilation cross section relative to the process XX ~ * 
A° — ► bb, which is presented in the left-hand panel of 
Fig. [3] The corrected cross section ctnlo is shown nor- 
malized to the leading order (LO) cross section clo, 
including either only the QCD diagrams, QCD and 
top-loop diagrams, or the full SUSY-QCD correction. 
The QCD correction is independent of the SUSY pa- 
rameters, while the top-quark loop correction is pro- 
portional to l/tan 2 /3 through the A°tt coupling. The 
points, for which the SUSY-QCD corrections are eval- 
uated, lie within the Higgs-funnel region and are in- 
dicated by a point in the left and centre panels of 
Fig. |U for /i < and fi > 0, respectively. It becomes 
clear that the annihilation cross section is decreased 
by about a factor two by the QCD contribution, where 
the principal contribution comes from the mass shift 
m,b — > rhb(s). The top-quark loop contribution, as al- 
ready mentioned, is negligible with respect to the QCD 
contribution at large tan/3 such as in the A -funnel re- 
gion. In the centre panel of Fig. [3] we see that A tov 
accounts for a few permille with respect to Z\qcd for 
tan/3 ~ 2, and only for less than 10 -5 for tan/3 > 40. 



When adding the SUSY-QCD contribution Z\susy, 
the annihilation cross section receives another impor- 
tant correction, see left-hand panel of Fig.[3j The SUSY 
correction decreases the cross section by another few 
percent at very low tan/3 and up to another 40 (10) 
percent for fj, < (fj, > 0), which underlines the im- 
portance of the full correction in the cosmological A- 
funnel region. The large difference between the correc- 
tion factors for fi > and /i < is mainly due to the 
mass renormalization, Eq. © , that is directly affected 
by sgn(/i). 

To evaluate the effect of the corrections on the cold 
dark matter relic density, we have implemented our full 
SUSY-QCD correction, as described above, into the 
DarkSUSY package [5] , which includes by default only 
the QCD corrections up to order C(a 2 ) without then- 
dependence on the energy scale s. Choosing either the 
pure Born or full SUSY-QCD calculation of the process 
XX ~~ * A, ~ > bb, we compare the resulting cold dark 
matter relic density i?CDM^ 2 to the observational 2a 
range in Eq. fTJ) to determine the favoured regions in 
the mo — mi/2 plane shown in the left-hand and centre 
panels of Fig. 2J As expected, the corrections do not 
affect the cosmological focus point (very low mi/ 2 ), 
bulk (low mo and mi/ 2 ) and coannihilation (low mo) 
regions, where the inspected process is suppressed with 
respect to other (co) annihilation channels. However, 
at large tan/3 = 44.5 (54) for fi < (n > 0), which 
still allow for electroweak symmetry breaking (EWSB) 
in a large region of the scanned mo — mi/2 plane, the 
LO allowed regions (light) are dramatically changed by 
the full SUSY corrections (dark), which reduce cr e ff as 
discussed above. The increase in i?cDM is compensated 
by a shift towards smaller SUSY masses. 

In the right-hand panel of Fig. [3] we plot the relic 
density for mo = 1200 GeV as function of m 1 / 2 - The 
graph shows the contributions of the different correc- 
tion terms. The effect of of the 0(a 2 s ) QCD corrections 
already included in DarkSUSY is considerably enhanced 
by our newly included 0(a 3 3 ) QCD and 0{a s ) SUSY- 
QCD corrections. Interesting effects due to the simul- 
taneous correction of the Higgs decay width Ta are ob- 
served, e.g. the minimum of J?cdm is shifted towards 
the resonance 2m x = itia- For a detailed discussion 
the reader is referred to Ref. [TT] . 

Finally we note a difference between the relic den- 
sity calculated with our code based on SPheno and 
FeynHiggs with respect to the one based on ISAJET, 
which is by default included in the DarkSUSY package. 
The right-hand panel of Fig. [4] shows the A-funnel re- 
gion for [i < 0, which is to be found at smaller values 
of mo. This difference should be due to the fact that 
the spectrum calculation is especially sensible to the 
spectrum generator for the neutralino and Higgs bo- 
son masses. In consequence, the resonance condition 
2m x = rriA, governing the A-funnel region, can be 
shifted in an important way. However, the effect of the 
SUSY-QCD corrections is the same, i.e. an important 
shift of the favoured region with respect to the leading 
order result. For more details on the impact of SUSY 
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Fig. 3. Left and centre: NLO annihilation cross section in units of the LO cross section, including the QCD, QCD and 
top, and the full SUSY-QCD correction and ratio of top-quark loop correction and QCD correction as function of tan f3. 
Right: Dark matter relic density based on LO, NLO QCD, and full NLO SUSY-QCD calculation as function of mi/2 
for mo = 1200 GeV, fi < 0, tan/3 = 44.5 (slope indicated in the left panel of Fig. HJ). 
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Fig. 4. Regions in the mSUGRA mo-mi/2 plane forbidden by a charged LSP / no EWSB and favoured by the observed 
cold dark matter relic density flcnwih 2 including LO (light) and full SUSY-QCD (dark) calculation. The corresponding 
SUSY spectrum has been obtained using SPheno / FeynHiggs for the left and centre panels, and using ISAJET for the 
right panel. 



spectrum calculations on dark matter annihilation see 
Ref. |22]. 



4 Conclusion 

We have presented a complete calculation of SUSY- 
QCD corrections to dark matter annihilation in the 
Higgs-funnel, resumming potentially large /^tan/3 and 
Aj, contributions. We have demonstrated numerically 
that these corrections strongly influence the extraction 
of SUSY parameters from cosmological data and must 
therefore be included in common analysis tools like 
DarkSUSY or micrOMEGAs. 
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